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EXECUTIVE  SUMMARY 


In  this  work,  we  present  a  method  to  estimate  the  acoustic  source  gain  of  targets  in  a  shallow-water 
environment,  using  passive  sonar.  This  is  accomplished  by  using  a  two-step  approach.  The  first  step  is  to 
localize  the  target,  which  is  accomplished  by  using  the  Sparsity-Cognizant  Source  Location  Mapping 
(scSLM)  algorithm,  developed  at  the  Space  and  Naval  Warfare  Systems  Center  Pacific.  This  algorithm  has 
been  shown  to  be  robust  to  environmental  mismatch  when  localizing  multiple  targets.  Moreover,  it  yields 
an  estimate  of  the  source  levels.  Unfortunately,  the  bias  in  these  estimates  can  be  large  depending  on  the 
values  of  the  algorithm's  tuning  parameters.  The  second  step  of  our  approach  is  introduced  to  decrease  the 
bias  in  the  acoustic  source  gain  estimates.  By  using  the  locations  obtained  by  scSLM  and  the 
environmental  model,  we  develop  an  estimator  that  yields  improved  acoustic  source  gain  estimates. 
Numerical  tests  arc  used  to  illustrate  the  quality  of  the  acoustic  source  gain  estimates  for  multiple 
broadband  targets  in  a  shallow-water  environment. 
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1.  INTRODUCTION 


Localization  of  an  underwater  target  and  then  estimation  of  its  acoustic  gain  using  passive  sonar  is  a 
difficult  task.  Localization  and  estimation  in  shallow-water  scenarios  is  particularly  challenging  due  to  the 
complexities  of  acoustic  propagation  through  this  environment.  For  low-frequency  signals,  acoustic 
propagation  in  shallow  water  is  characterized  by  multipath  due  to  the  multiple  interactions  that  they  sustain 
with  the  sea  surface  and  sea  floor,  by  refraction  due  to  a  depth-dependent  sound  speed  profile,  and  by  the 
spatio-temporal  variability  of  environmental  parameters.  Shipping  lanes  and  marine  life  also  contribute 
noise  that  masks/interferes  with  the  acoustic  signals  of  interest  and  leads  to  low  signal-to-noise  ratios 
(SNRs)  at  the  hydrophone  array  [1], 

Various  matched-field  processing  (MFP)  algorithms  are  useful  in  localization  of  underwater  acoustic 
sources.  MFP  assumes  that  sources  arc  located  in  a  predefined  area  and  uses  an  acoustic  propagation  model 
to  predict  the  pressure  fields,  known  as  replicas,  at  the  hydrophone  array  for  a  grid  of  locations.  These 
replicas  arc  then  ’’matched”  to  the  hydrophone  measurements  to  produce  an  ambiguity  surface,  which 
summarizes  the  acoustic  power  estimates  across  the  grid.  Then  localization  is  achieved  by  choosing  the 
largest  peaks  on  the  ambiguity  surface  and  assigning  them  to  source  locations.  Classical  MFP  suffers  from 
large  sidelobes,  which  makes  localization  of  multiple  sources  difficult  and  can  be  sensitive  to 
environmental  mismatch  [2]. 

Different  from  MFP,  Sparsity-Cognizant  Source  Location  Mapping  (scSLM)  casts  the  localization 
problem  as  a  group-Lasso  regression  problem.  It  has  been  shown  that  location  estimates  obtained  via 
scSLM  arc  more  accurate  than  the  ones  obtained  by  classical  MFP  algorithms,  such  as  Bartlett  and  Capon 
[3].  It  has  also  been  shown  that  they  arc  robust  to  environmental  mismatch  [3].  This  work  uses  the  scSLM 
algorithm  to  localize  the  targets  and  estimate  their  acoustic  gain.  It  is  then  shown  that  the  acoustic  source 
gain  estimates  provided  by  the  scSLM’s  output  do  not  accurately  estimate  the  true  source  gain.  A  two-step 
approach,  using  scSLM  for  the  localization  and  a  least-squares  estimator  for  obtaining  the  source  gains 
yields  improved  source  gain  estimates.  This  method  differs  from  ones  found  in  literature,  that  tend  to  use  a 
simplified  passive  wave  equation  to  estimate  the  acoustic  source-levels  of  ships  and  marine  life  [4-6].  One 
of  the  main  problems  with  using  a  simplified  passive  wave  equation  in  a  shallow-water  environment  is  the 
fact  that  it  does  not  account  for  multipath  or  spatio-temporal  variability  of  environmental  parameters. 

Notation  :  upper  (lower)  boldface  letters  are  used  for  matrices  (column  vectors);(-)’ ((•)')  denotes 
matrix  and  vector  transpose  (conjugate  transpose); [-]n,m([']n)  the  (n,m)-entry  of  a  matrix  (n-entry  of  a 

vector);  ||y||?  =  (E^r  \Vn\q)  ^ 
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2.  PROBLEM  FORMULATION 


An  array  of  N  hydrophones  with  known  and  arbitrary  geometry  is  used  to  capture  T  discrete-time 
acoustic  measurements  {y(t)}J=1  from  K  sources.  Let  yn(t)  :=  \y(t)\„  coiTespond  to  the  acoustic 
measurements  gathered  by  the  n-th  hydrophone  at  time  t.  With  hn(t:  r/.)  being  the  discrete-time  impulse 
response  of  the  channel  at  time  t  and  unknown  source  location  17.  G  M.d,  where  d  G  {2,  3},  yn ( t )  is 
modeled  as 

K 

Unit)  =  y sk )  ( t ;  rfc)  +  en(t ),  n  =  l,...,N  (1) 

fc= 1 

where  >1=  denotes  the  convolution  operator,  .sy.  (t)  the  acoustic  signature  for  the  k-th  source  and  en (t)  is  a 
zero-mean  random  measurement  noise. 

Although  it  is  possible  to  work  directly  in  the  time  domain,  the  computational  complexity  required  to  do 
so,  can  be  excessive  [7].  Instead,  the  problem  is  cast  in  the  frequency  domain.  First,  {y(t)}J=1  is 
partitioned  into  M  blocks  and  transformed  to  the  frequency  domain  via  the  discrete  Fourier  transform 
(DFT).  Per  frequency  w j,  let  ym,f  G  CN  and  s^j  G  C  denote  the  DFT  coefficient  vector  for  the  m-th 
measurement  block  and  the  acoustic  gain  DFT  coefficient  for  the  n-th  source,  respectively.  Then, 

{y m,f}Z=\  is  modeled  as 

K 

y rn,f  =  y 1  Sk,f  (p k,f  +  Vfc,/)  +  em,f ,  m  =  1, ...,  M  (2) 

fc= 1 

where  emj  denotes  the  Fourier  coefficients  at  Wf  corresponding  to  the  noise  in  the  m-th  block;  p^.  jeCN 
the  Wf  model-predicted  Fourier  coefficients  at  the  array  for  a  source  located  at  ry.,  which  will  be  referred  to 
as  replicas;  and  v/,.j  G  CN  an  unknown  perturbation  vector  effecting  Pkj,  which  helps  account  for  the 
mismatch  between  the  acoustic  propagation  model  and  the  true  propagation  environment. 

Let  Y f  :=  [yqj, ....  Ymj]  comprise  all  M  Fourier  coefficients  vectors  at  Wf.  Given  an  underwater 
propagation  model,  Fourier  coefficients  {Y  f}f=\  for  F  frequencies  {  w  j  }  ^  f  and  the  number  of  sources 
K.  The  localization  process  seeks  to  find  the  positions  {rj.}*=1  while  still  being  robust  to  unknown 
perturbations  {{vfcj}^=1}J_1.  Even  if  the  adopted  underwater  propagation  model  were  to  exactly 
characterize  the  propagation  environment,  solving  for  {rfc}^=1  is  difficult  because  it  entails  a  nonlinear 
regression  problem.  A  closed-form  expression  relating  17.  to  p/.j  is  unavailable  since  in  most  cases  of 
interest,  finding  p/.j  entails  solving  the  wave  equation  for  specific  boundary  conditions  given  by  the 
environment  model  [1], 

2.1  LOCALIZATION 

The  scSLM  algorithm  circumvents  some  of  the  computational  challenges  posed  by  the  nonlinear 
relationship  between  ry.  and  {p/,,/}|=P  by  constructing  a  grid  Q  :=  { ry,  G  that  spans  the  region  of 

interest.  Let  p gj  denote  the  replica  vector  at  frequency  Wf  for  a  source  located  at  rg.  For  a  given  Wf, 

P f  :=  [pi.f , ....  Pa./}  G  C;Vx6'  comprises  all  normalized  replicas  corresponding  to  Q  where 
||Pg,/ II2  =  L  Vry.  Each  y m,f  is  now  modeled  as 

G 

y  rnj  =  S9,f(PgJ  +  v9,/)  +  em,f’  /  (3) 

9=1 

where  v9  j  G  C^'  denotes  the  unknown  perturbation  vector  affecting  p gj.  Since  K<G,  most  of  the 
{sgj}’s  are  expected  to  be  zero. 
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Let  s  f  :=  [SlJ, ...,  sgj}'  £  CG  denote  the  complex-valued  vector  of  regression  coefficients  at  Wf, 

S  :=  [si ,  .4;  sp]  £  CGxF  denote  the  complex-valued  matrix  comprising  all  regression  coefficients, 

V/  :=  K /,-  ,  V(;,r]  £  CNxG  a  complex- valued  matrix  comprising  all  perturbation  vectors.  The  scSLM 
algorithm  uses  an  t\/lq  regularizer,  also  known  as  a  group-Lasso  regularize!',  to  find  the  optimal  S  while 
being  robust  to  environmental  mismatch.  Equation  3  can  be  cast  as  a  convex  optimization  problem  that 
forces  sparsity  and  coherence,  in  the  source  locations  across  frequencies,  see  [8,  9]  for  a  detailed 
description.  In  this  work  coherence  refers  to  a  common  support  across  source  localization  maps  at  different 
frequencies,  which  is  reasonable  to  exploit  since  it  is  assumed  the  source  locations  arc  immutable  across 
frequencies.  The  scSLM  algorithm  solves  the  following  minimization  problem: 


1 


F 


min  -  >  E 

secGxF  2 M  ^ 


Yf-(Pf  +  Vf)sflM 


+  A4 


G 

E 

0=1 


(4) 


where  qg  is  defined  as  the  g-th  row  vector  of  S,  qg  :=  [s9,i, ...,  s9,f]-  //  >  0  is  a  tuning  parameter 
controlling  the  sparsity  level  on  S  and  q  £  (1, 00].  The  scSLM  algorithm  assumes  that  the  columns  of  Vy 
arc  zero-mean  random  vectors  with  covariance  matrix  £  /■  <G  CGxG  and  independent  across  V  f  s. 

By  choosing  //,  one  can  construct  a  map  with  K  nonzero  locations  that  correspond  to  the  source 
locations  [10]. 


2.2  SOURCE  LEVEL  ESTIMATION 


One  should  note  that  Equation  (4)  gives  both  an  estimate  of  the  source  location  and  the  acoustic  gains. 
However,  since  sgj  was  calculated  using  normalized  replicas  the  estimate  does  not  correspond  to  the  true 
acoustic  gain,  but  rather  it  corresponds  to  the  acoustic  source  gain  decreased  by  the  acoustic  propagation 
loss  from  the  source  to  array.  To  compute  the  source  gain  estimate,  we  need  to  remove  the  normalization  of 
the  replica  that  was  used  to  obtain  sgj.  Letting  {()l  )G  \  be  the  grid  points  that  correspond  to  the  K  top 
ranked  intensities  from  the  scSLM  algorithm,  we  obtain  the  following  acoustic  gain  estimates 

S9lJ 

S9l,f  ~ 

P  9t,f 


(5) 


One  caveat  with  this  method  of  estimating  the  source  level,  sg*j  is  that  it  yields  biased  estimates  since  it  is 
based  on  the  group-Lasso  regularizer  [8]. 


To  circumvent  this  problem,  we  propose  to  use  the  estimated  location  of  the  source  g*k,  but  will  then  use 
the  least-squares  estimator  to  find  the  source  gains.  Since  the  scSLM  algorithm  guarantees  a  common 
support  and  each  frequency  has  an  independent  source  level,  we  can  obtain  the  source  gains  for  each 
frequency  separately.  Letting  Pj>  :=  [p9*,/, ...,  Pg*K,f]  and  sF  :=  [s9*,/,  •••,  sg*Kj]',  we  postulate  the 
following  estimator  for  the  source  levels,  sj: 


~K  ■  II  t>K  K 

si  :=  argrnm  yy  —  P  f  s  f 

sKecK><  1 

which  has  the  following  closed  form  solution: 


2 

2  ’ 


s/ 


-1 


y  /• 


(6) 

(7) 


Note  that  if  there  is  only  one  source  present.  Equation  (7)  simplifies  to  the  following  form: 


s/ 


p  g1*,/ty/ 

P 9i’f  2 


(8) 
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3.  NUMERICAL  TESTS 


This  section  serves  to  illustrate  the  estimation  quality  of  the  acoustic  source  gains  in  a  broadband 
multi-source  environment,  using  passive  sonar.  Numerical  tests  on  synthetic  data  sets  were  conducted  to 
illustrate  the  performance  of  the  proposed  algorithm. 

A  shallow- water  environment,  as  depicted  in  Figure  la,  was  created.  This  environment  contains  a 
vertical  line  hydrophone  array,  which  consists  of  N  =  25  elements  that  span  the  entire  water  column  of 
100  m.  A  grid  of  10,25 1  locations  spanning  a  range  0  to  10  km  and  depths  0  to  100  m  was  used.  This  gives 
a  range  and  depth  spacing  of  50  and  2  m,  respectively.  All  replicas  for  this  environment  were  calculated 
using  KRAKEN  [11]. 

We  first  consider  a  single  broadband  source  and  the  following  five  frequencies:  65,  150,  220,  300,  and 
350  Hz.  These  five  frequencies  correspond  to  the  following  source  gains:  95,  97,  100,  97,  and  95  dB, 
respectively  (Figure  lb).  We  then  introduce  additive  white  Gaussian  noise  to  the  replicas  to  obtain  a 
predefined  average  SNR  across  the  frequencies.  The  average  SNR  was  calculated  at  the  closest  point  of 
approach  (CPA)  and  will  vary  according  to  the  targets  distance  (Figure  2). 

Next  we  simulate  a  two  source  case.  Both  sources  are  emitting  at  the  same  five  frequencies  and  same 
acoustic  gains:  65,  150,  220,  300,  and  350  Hz,  and  95,  97,  100,  97,  and  95  dB,  respectively.  The  difference 
in  the  two  sources  arc  the  CPA  and  depth  values  as  can  be  seen  in  Figure  3. 

For  all  tests,  we  considered  a  correct  localization  to  be  a  square  of  ±1  around  the  true  location,  as 
shown  in  Figure  4.  We  also  do  not  allow  multiple  targets  to  localize  within  the  same  area  of  localization. 
For  the  following  tests,  scSLM  with  q  =  2  was  used. 

3.1  RESULTS 

The  following  results  show  the  average  estimates  over  50  Monte  Carlo  runs.  scSLM  parameters  a  and  A 
were  set  to  0.005  and  0.85  respectively. 

3.1.1  ONE  SOURCE 

For  the  single  source  case,  the  average  SNR  at  CPA  was  set  to  3  dB.  The  resulting  SNR's  per  frequency 
is  shown  in  Figure  2.  Figure  6  shows  that  the  scSLM  acoustic  source  gain  estimations  are  not  an  accurate 
portrayal  of  the  true  signal  strengths.  Again,  this  is  due  to  the  innate  bias  of  scSLM,  that  stems  from  its 
reliance  on  the  group-Lasso  regularizer.  Since  this  bias  occurs  in  all  the  scSLM  source  level  estimations, 
we  will  no  longer  show  these  results.  When  first  localizing  the  sources  based  on  scSLM  and  then  using  the 
least-squares  estimator  from  Equation  (6)  to  estimate  the  source  levels,  we  arc  able  to  better  approximate 
the  true  acoustic  gain  (Figure  7).  Note  that  as  the  SNR  drops,  localization  becomes  more  difficult.  Hence, 
estimating  the  source  acoustic  gains  also  becomes  more  difficult.  This  is  shown  in  the  larger  error-bars  in 
Figures  6  and  7. 
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(b)  Signal  power 


(c)  True  source  -  Range 


(d)  True  source  -  Depth 


Figure  1.  Simulation  setup  -  Schematic  of  the  underwater  acoustic  propagation  model  used  for  the 
simulation.  Sample  parameter  values  for  this  acoustic  propagation  model  are  14  =  1, 500  m/s,  v2  =  1, 460 

m/s,  v3  =  1, 600  m/s,  v4  =  1, 750  m/s,  v5  =  1, 750  m/s,  ab  =  0.2  dB/m/KHz,  and  pb  =  1.8  g/cm3. 


Time 


Figure  2.  SNR  of  each  frequency  -  one  source. 
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(a)  True  source  -  Range  (b)  True  source  -  Depth 

Figure  3.  Simulated  range  and  depth  of  the  two  sources. 


Figure  4.  Correct  localization  region. 


(a)  Range  (b)  Depth 

Figure  5.  Location  estimates  from  scSLM  for  a  single  source. 
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(a)  Frequency  =  65  Hz 
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(c)  Frequency  =  220  Hz 


(d)  Frequency  =  300  Hz 


Time 


(e)  Frequency  =  350  Hz 


(f)  Average  error 


Figure  6.  Single  Source:  Acoustic  source  gain  estimates  obtained  via  scSLM.  The  stem  plots  represent 
the  percentage  of  correct  localizations  estimated  via  50  Monte  Carlo  runs.  The  source  gain  plot  shows  the 
average  gain  magnitude  in  dB,  20  logw(\  sj),  along  with  its  standard  deviation.  The  average  error  plot 
corresponds  to  ef  =  20  logw(\\sf\  -  |  SJ|).  The  true  source  gain,  at  frequency  /,  and  the  average  source 
gain  from  the  Monte  Carlo  runs  is  s/,s/,  respectively. 
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(a)  Frequency  =  65  Hz  (b)  Frequency  =  150  Hz 


(c)  Frequency  =  220  Hz 


(d)  Frequency  =  300  Hz 


(e)  Frequency  =  350  Hz 


(f)  Average  error 


Figure  7.  Single  Source:  Acoustic  source  gain  estimates  obtained  via  least-squares.  The  stem  plots 
represent  the  percentage  of  correct  localizations  estimated  via  50  Monte  Carlo  runs.  The  source  gain  plot 
shows  the  average  gain  magnitude  in  dB,  20logj0(\  5J),  along  with  its  standard  deviation.  The  average 
error  plot  corresponds  to  ef  =  20  fogi0(||s/|  -  |  sj  |).  The  true  source  gain,  at  frequency  /,  and  the  average 
source  gain  from  the  Monte  Carlo  runs  is  sf,sf,  respectively. 


Figure  8.  Location  estimates  for  two  sources. 


Figure  9.  SNR  for  each  target. 


3.1.2  TWO  SOURCES 

In  the  two-source  simulation,  we  increased  the  average  SNR  at  CPA  to  8  dB  (Figure  9).  Target  1  will  be 
defined  as  the  source  at  the  depth  of  52  m,  while  target  2  will  be  the  source  at  72  m  (Figure  3).  The  results 
of  the  two-source  Monte  Carlo  simulation  are  shown  in  Figures  10  and  11.  These  figures  illustrate  that 
correct  localizations  on  multiple  targets  can  be  difficult,  even  at  higher  SNRs  (compared  with  the 
single-source  case).  These  figures  also  show  that  target  1  was  localized  more  often  than  target  2  and  that 
the  error-bars  on  target  1  at  each  frequency  do  not  vary  as  much  as  target  2.  The  smaller  variations  are 
likely  caused  by  the  incorrect  localization  providing  fallacious  source  gain  estimates.  As  the  targets 
approach  CPA,  the  acoustic  gains  were  closer  to  the  expected  values,  along  with  having  lower  variances. 
This  is  most  likely  due  to  the  increase  in  localization  percentages  and  higher  SNRs. 
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(a)  Frequency  =  65  Hz  (b)  Frequency  =  150  Hz 


(c)  Frequency  =  220  Hz 


(d)  Frequency  =  300  Hz 


Time 


(e)  Frequency  =  350  Hz 


(f)  Average  error 


Figure  10.  Two  Sources:  Target  1  acoustic  source  gain  estimates  obtained  via  least-squares.  The  stem 
plots  represent  the  percentage  of  correct  localizations  estimated  via  50  Monte  Carlo  runs.  The  source  gain 
plot  shows  the  average  gain  magnitude  in  dB,  20  logw(\  sj),  along  with  its  standard  deviation.  The  average 
error  plot  corresponds  to  ef  =  20  log10(\\sf\  -  |  sj|).  The  true  source  gain,  at  frequency  /,  and  the  average 
source  gain  from  the  Monte  Carlo  runs  is  sf,sf,  respectively. 
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(e)  Frequency  =  350  Hz 
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Figure  11.  Two  Sources:  Target  2  acoustic  source  gain  estimates  obtained  via  least-squares.  The  stem 
plots  represent  the  percentage  of  correct  localizations  estimated  via  50  Monte  Carlo  runs.  The  source  gain 
plot  shows  the  average  gain  magnitude  in  dB,  20  logw(\  sj),  along  with  its  standard  deviation.  The  average 
error  plot  corresponds  to  ef  =  20  log10(\\sf\  -  |  sj|).  The  true  source  gain,  at  frequency  /,  and  the  average 
source  gain  from  the  Monte  Carlo  runs  is  sf,sf,  respectively. 
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4.  CONCLUSIONS 


This  work  focused  on  estimating  the  acoustic  source  gains  of  broadband  sources  in  a  shallow- water 
environment.  The  results  show  that  accurate  acoustic  source  gain  estimation  is  possible  if  the  sources  were 
correctly  localized.  Knowledge  about  the  environment  was  a  fundamental  enabler  for  the  acoustic  source 
gain  estimation.  When  dealing  with  multiple  sources  the  source  gain  estimation  does  not  perform  as  well, 
but  if  localization  was  achieved,  then  it  is  still  able  to  give  accurate  results  at  higher  SNRs.  Future  research 
will  explore  how  this  algorithm  performs  on  real  data  set(s)  and  in  the  presence  of  environmental 
uncertainty. 
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